IDLA ON THE SUPERCRITICAL PERCOLATION CLUSTER 



ERIC SHELLEF 

Abstract. We consider the internal diffusion limited aggregation (IDLA) 
process on the infinite cluster in supercritical Bernoulli bond percolation on 
7L d . It is shown that the process on the cluster behaves like it does on the 
Euclidean lattice, in that the aggregate covers all the vertices in a Euclidean 
ball around the origin, such that the ratio of vertices in this ball to the total 
number of particles sent out approaches one almost surely. 



1. Introduction 

1.1. Background and discussion. Given a graph, IDLA defines a random ag- 
gregation process, starting with a single vertex and growing by a vertex in each 
time step. To begin the process, we specialize a vertex v to be the initial aggregate 
on the graph. In each time step, we send out one random walk from v. Once this 
walk exits the aggregate, it stops, and the new vertex is added to the aggregate. 
Let I(n) be the aggregate in step n, thus 1(1) = {v} and = n. 

This process is a special case of a model Diaconis and Fulton introduced in 
[DF91]. In the setting where the graph is the d-dimensional lattice, Lawler, Bram- 
son and Griffeath in |LBG92| used Green function estimates to prove the process 
has a Euclidean ball limiting shape. Let B r be all vertices in 7L d of Euclidean 
distance less than r from the origin. The main theorem in [LBG92 implies that 
for any e > 0, I(\Bn\) D £>(i_ c )a for all sufficiently large R with probability one. 
Seeking to generalize this result to other graphs, we note that convergence of a 
random walk on the lattice to isotropic Brownian motion plays an important role 
in convergence of IDLA on the lattice to an isotropic limiting shape. 

However, this property by itself is in general not enough for IDLA to have such 
an inner bound, as the following example shows. Consider the three dimensional 
Euclidean lattice, and choose a vertex v at distance R from the origin. Let r = R 9 
and let B r (v) be the set of vertices of distance less than r from v. Remove all edges 
but one from the boundary of B r (v), and denote by v' the only vertex in B r (v) with 
a neighbor outside of B r (v). Let us look at X(\B2r\). A rough calculation gives that 
the average number of visits to v' is of order |-E?2ii|- Since at least |_B r (w)| = R 27 
visits to v' are needed to fill B r (v), we don't expect the ball to be full after an order 
of R 3 particles have been sent out. Repeating this edge removal procedure for balls 
of radius R% 9 at distance R n from the origin where R n — 2 n R, will ensure that 
there is never a Euclidean inner bound. However, a random walk from the origin on 
this graph will converge to Brownian motion because our disruptions are sublinear. 
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FIGURE 1.1. IDLA of 1000 particles from marked vertex on the per- 
colation cluster with p — 0.7. 



We will not give full proofs of these facts, but hope they convince the reader that 
to get an inner bound some kind of local regularity property is needed. 

The main theorem of this paper states that IDLA on the supercritical cluster 
in the lattice has a Euclidean ball as an inner bound. The two main tools that 
are used to show this are a quenched invariance principle [BB07 which gives us 
convergence in distribution to Brownian Motion from a fixed point, and a Harnack 
inequality from [Bar04 , which give us oscillation bounds on harmonic functions 
in all small balls close to the origin. The latter allows us to establish the local 
regularity missing from the above example. 

1.2. Assumptions and statement. Consider supercritical bond percolation in 
Z d with the origin conditioned to be in the infinite cluster. Let the graph r(V, E) 
be the natural embedding in M. d of the infinite cluster, i.e. V C M. d . Fixing this 
embedding for T, we get two separate (but comparable, see |AP96j ) distances. We 
denote by \x — y\ or d(x, y) Euclidean distance between points x, y, and by dr{x, y) 
the graph distance between them. If one of the points is not in V, dr(x,y) = oo. 
Let B r (x) = {v € V | d(x, v) < r} be the vertices contained in a Euclidean ball of 
radius r and center x. We abbreviate B r (0) as B r . To differentiate between such 
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a set of vertices and a full ball in R d , we denote by bold lowercase the following: 
h r (x) = {y E R d | d(x,y) < r}. Let d r (x) be a box of side r with center x, and as 
above, let D r (x) be the vertices in this box. 

By existing work, which we will reference later, we know that with probability 
one, the graph of the supercritical cluster, T, in its natural K d embedding, satisfies 
the following assumptions: 

(1) A random walk on T converges weakly in distribution to a Brownian Motion 
as defined in II .41 



(2) Convergence to vertex density ag > as defined in 1.5 

(3) A uniform upper bound on exit time from a ball as defined in 1.6 

(4) A Harnack inequality as defined in |1.7| 

We denote by X t a discrete "blind" random walk on V defined as follows. For 
x G T and ( e NU {0}, P (X t+1 = y\X t = x) = l/2d if y is a neighbor of x in 
F, and P {X t +\ = x\X t = x) = 1 — degx/2d. We prefer this walk since its Green 
functions are symmetric, a fact which will be useful later on. For non-integer t we 
set X t = X\t\ ■ 

Note that assumption [4] does not imply [3] see |Del02| for a counterexample. 

Using only above assumptions and that V C Z d (which serves mostly to simplify 
notation and could be replaced by weaker conditions), we show the IDLA process 
starting at will have a Euclidean ball inner bound as stated in the following 
theorem: 

Let I(n) denote the random IDLA aggregate of n particles starting at 0. 

Theorem 1.1. Almost surely, for any e > 0, we have that for all large enough R, 
B {1 _ e)R cX(\B R \) 

We fix e with which we prove the above for the rest of the paper. 

1.3. Outline. In the remainder of this section we state our assumptions precisely, 
and explain why these assumptions are valid almost surely for the infinite cluster 
in supercritical percolation. 

Then, to prove IDLA inner bound, we need to show that for each vertex v G V 
the Green function from and expected exit time from a ball around of a random 
walk starting at u, behave similarly to those functions of a Brownian motion. The 
exact statement needed appears in Lemma |3.1| The invariance principle gives us 
integral convergence of these functions to the right value, but to improve them 
to pointwise statements, we must show that they are locally regular. We use the 
Harnack inequality from |Bar04| to prove they are Holder continuous. This is a 
similar scheme to improving a CLT to an LCLT, as was done in [BH09 using a 
different method than ours. 

In section [2] we start by proving a lemma comparing expected exit time from a 
set to the Green function of a point in the set. Then, we show how our Harnack 
assumption leads to an oscillation inequality which we use to show regularity of the 
Green function and expected exit time of the walk in a ball when we are far from the 
boundary and the center. Next, in section|3] we use our assumption of an invariance 
principle to show that in small balls the integral of these functions approaches a 
tractable limit that can be calculated by knowledge of Brownian motion. Finally, 
in section [4] we utilize these estimates to prove theorem |1.1| 

Remark 1.2. Since the paper makes use of results based on ergodic theory, there is 
no rate estimate for the convergence in theorem |1.1| Also, on the Euclidean lattice 
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there is an almost sure outer bound with sublinear fluctuations from the sphere. 
Another interesting question that is not treated here, is whether a similar outer 
bound holds in our setting. 

1.4. Weak convergence to Brownian Motion. In this and the next three sub- 
sections, we give a precise formulation of each of our assumptions from above, and 
argue that they hold a.s. on T, the supercritical cluster. 

Let = C([0,T] — > Mr), i.e. the continuous functions from the closed interval 
[0, T] to (i-dimensional Euclidean space. For R € N, let 

w R (t) = ~ (X [tR2] + (tR 2 - [tR 2 \)(X [tR 2 ]+1 - X vm ^)) . 



Thus Wfi(t) is a scaled linear interpolation of X t (defined in 1.2 1 with its restriction 
to [0,T] an element of C%. 

We say that assumption [I] holds if for any T > 0, the law of w R (t) on con- 
verges weakly in the supremum topology to the law of a (not necessarily standard) 
Brownian Motion {B t : < t < T) . 

Lemma 1.3. Assumption r7] holds for T with probability one. 

Proof. This is Theorem 1.1 in [BB07 . Another paper with a similar invariance 
principle is |MP05| . □ 

Let B(t) denote the Brownian motion weak limit of X t . 

1.5. Convergence to positive density. Assumption [2] holds if there exists a 
positive ac such that for any 5, 7 > and all sufficiently large R we have that for 
all x e b R and all SR < r < R 

\B r {x) 



\Dr(x)\ 

a G 



< 7- 



That is, balls and boxes of size 5R to R, have vertex density in (ag — 7, og + 7) 
for all large enough R. 

Lemma 1.4. Assumption [1| holds for T with probability one. 

Proof. Let 9(p) be the probability for to be in the infinite cluster. 8(p) is positive 
in the supercritical regime. From Theorem 3 in [DS88 for d = 2 and from (2) on 
p. 15 of [Gan89 for d > 2, we know that for any p > there is a positive c = c(p) 
such that 

P (n- d \D n (x)\ < (1 - p)9(p)) < exp(-cn). 

where x g M. d and n E N. Recall D n (x) are the vertices in a box of side-length n 
and center x. Theorem 2 in |DS88| proves the easier density upper bound for d = 2, 
which under trivial modifications gives that for all d 

P (n~ d \D n {x)\ > (1 + p)0{p)) < exp(-cn). 

The above, together with Borel Cantelli gives that if we choose 7 > 0, x G bi, then 
for all large enough R 

\\D pR {Rx)\-6{p){pR) d \ < ^R d . 

So we have the result for small boxs of size pR. To expand it to larger boxes, let D 
be any box of diameter between 8R to R. We partition the box D R into pi?-sized 
boxes and choose p = p(S) small enough so that the number of boxes that intersect 
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both D and Dr\D is negligible compared to the number that intersect D. Proving 
for balls is similar. □ 

1.6. Uniform Bound on Exit time from ball. Denote by T r (x) the first time 
the walk leaves B r (x), and write r r for r r (0). 

Assumption [3] holds if there is a ce such that for any 6 > and all large enough 
R, if x E Br and r > 5R then 

(1.1) E x [r r (x)}<c E r 2 . 

The assumption holds for V with probability one as a direct consequence of the 
following lemma, proved below. 

Lemma 1.5. With probability one, there is an L such that for all 5 > and all 
large enough R, if x G Br and r > SR then 

(1.2) P x (T r {x) > Lr 2 ) < exp(-cL). 

Proof. Since we prove for all R outside a bounded interval, it suffices to prove the 
above for r = SR with fixed S > 0. To prove this holds a.s. for supercritical 
percolation clusters, we use a heat kernel upper bound given in Theorem 1 of 
Barlow's paper [Bar04| . The proof in |Bar04| is for continuous walks with a mean 
time of one between jumps. However, it can be transferred to our discrete walk 
using Theorem 2.1 of [BB07 . We state an implication of what was proved. 

With probability one there exists a function from V to N, \l~ x < °°} x< zv where 
T x is sublinear in the sense that tSt — > 0, and there exist positive constants c\,c 2 
such that for any x E V, the following holds. For all y E B t (x), t > T x - 

(1.3) P x (X t =y)< c 2 t- d ' 2 . 

Barlow's result actually states that almost surely T x grows slower than a logarithmic 
function of ||x|| and gives stronger Gaussian bounds for the heat kernel from below 
and above. 

Using the heat kernel upper bound (1.3), the convergence to zero of t?%, and 



the upper bound on vertex density resulting from T C Z d , we have that for any K 
and all large enough R, if x G Br and r = 5R, then 

P x {X Kr .2 £ B r {x)) = 1- J2 Px{X Kr .2=y) 

yeB r (x) 

> 1- Y, c 2 K- d ' 2 r~ d 

yGB r (x) 

> 1 - C(d)K- d/2 . 

Thus for some K' , for all large enough R, and all x E Br, P x {Xk' T 2 G B r (x)) < |. 
Next, for any positive L, we use the Markov property to upper bound P x (r r {x) > 
LK'r 2 ) by 

VL\ 

P x (X K , r 2 E B r (x)) ]JP(X iK , r3 E B r (x) | X (i _ 1)K , r 2 E B r (x)) < 

i=2 

which proves the claim. □ 

We will later use that assumption [3] implies L\ convergence of R~ 2 tr A T to 
R~ 2 tr with T, uniformly in R for all x E Br. 
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Lemma 1.6. For all f3 > there is a T{j3) such that for all large enough R, 
Vx e B R 



(1.4) 



\E x [t r ]-E x [t r ATR 2 ]\ <f3R 2 



We rewrite the left hand side of (1.4) and use the Markov property with the exit 
time assumption. 



E x [(t r -TR 2 )1 {tr>tr2} ] = 

y£B R 



E x [t r - TR 2 | Tfl > TR 2 ,X{TR 2 ) = y] 



P(T R >TR 2 ,X(TR 2 )=y) 
= E v t r «] > TR 2 ,X{TR 2 ) = y) 

< c E R 2 P ( T R > TR 2 ,X{TR 2 ) = y) 
y£B R 

= c e R 2 P{t r >TR 2 ). 

Since E x [t r /R 2 ] is bounded by ce for all R, we have by the Markov inequality 
that P(t r /R 2 > T) < c E /T which converges to zero as T goes to oo uniformly in 
R for all x £ B R . 



1.7. Harnack inequality. A function h : T — > M is harmonic on x E T if h(x) = 
l j ^( a; )r 1 J2 y eN(x) h(y) where N(x) = {y G T : d r (x,y) = 1} are the neighbors of 
x. We say h is harmonic onZcT if it harmonic on every vertex v G Z. 

Assumption [4] holds if there is a 1 < ch < 00 such that for all 5 > we have that 
for all large enough R, if x € B Rl r > SR and h is a function that is non-negative 
and harmonic on B2 r {x) then 

(1.5) sup h < ch inf h. 

B r (x) Br{?) 

In order to simplify the paper, the inequality above is formulated for Euclidean 
distance rather than graph distance. We show that since the distances are almost 
surely comparable on T, this is the same. We start by proving the assumption 
holds for graph distance, using Lemma 2.19 and Theorem 5.11 of [Bar 04 . We 
write B^ = {x £ T : dr(0,x) < r)}. The lemma, using appropriate parameters and 
Borel Cantelli, tells us that for all large enough R, B RAogR are very good balls. 
Specifically, all balls contained in B^ loa . R of graph radius larger than _R 1 / 4 ( d + 2 ) 
have a positive volume density and satisfy a weak Poincare inequality as explained 
in (1.15) and (1.16) of the same. Theorem 5.11 then tells us that for very good 
ball of graph radius RlogR, all x £ ^(ijiog_R)/3 sa tisfy that for any function h 
non-negative harmonic on B R (x) 

(1.6) sup h < ch inf h 

where CH(d,p) > is a constant dependent only on dimension and percolation 
probability. Since all but a finite number of balls of graph radius R log R are very 
good and R is o(i? log R) we have assumption [4] for graph distance. 
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Next, we transfer this to the Euclidean balls formulation of (1.5). By Theorem 



1.1 of |AP96| . we have that for some k > 0, p(d,p) > and M < oo 

dr{x, y) > pm x, y € T, \x — y\ = m > M < exp(— km). 

Let A(x,y) be the event {dr{x,y) > pm}. Union bounding the probability for 
A(x,y) over every pair of points in Br of Euclidean distance greater than clogi? 
for some large c(fc), we upper bound the probability of any such event occuring by 



CR d md_1 exp(-fcm) < C ' R d {\ogR) d ^ RT 2a 

m—c log R 

which is summable for d > 1. Hence by Borel Cantelli, almost surely for all large 
enough R, we have that for any x,y € Br with \x — y\ > clogi?, dr(x,y) < 
p\x — y\. Since logi? is o(R), this implies that for any 5 > and all large enough 
R, for any x € Br, 

(1.7) B Sp -t R (x) c B t S r{x). 

Note that B^(x) C B r (x) always because Y is embedded in Z d . Hence, given a 
function that is non-negative harmonic on B2 r (x) it is also such on B\ r (x). For 
large enough R, we use (1.6| and (1.7) to get (1.5) for B p -i r (x) with constant ch- 



A routine chaining argument (see e.g. (3.5) in |Del02j) transfers this to B r (x) as 
required with new constant c#. 



2. POINTWISE BOUNDS ON GREEN FUNCTION AND EXPECTED EXIT TIME 

In this section, we show the assumptions of a uniform bound on exit time from 
a ball along with the Harnack inequality give us pointwise bounds of the Green 
function and expected exit time of a random walk. 

As a convention, we use plain c and C to denote positive constants that do 
not retain their values from one expression to another, as opposed to subscripted 
constants Cj, that do. In general, these constants are graph dependent, and in the 
context of percolation, can be seen as random functions of the percolation configu- 
ration. However, we view the graph as being fixed and satisfying the assumptions 
stated in subsection 11.21 

We start with a general lemma on the relation between expected exit time from 
a set and the expected number of visits to a fixed point in the set. 

Let Z C r and let tz be the first hitting time of Z for X t . For x € Y we set 
G z {x) = E x o l{A- t =x}J i the expected number of visits to x of a walk starting 
at x before tz- 

Lemma 2.1. There is a k — k(d) > such that for any x € Y and Z C Y where 
N(x) n Z = 



(2.1) 



E x [Tz]>kG 2 z /\ogG z 



Proof. Recall from |1. 2 1 that X t has a positive staying probability at certain vertices. 
Since we apply electrical network interpretation to estimate hitting probabilities, 
we prove the above for Y t , the usual discrete simple random walk, with probability 
to stay at a vertex. This implies (2.1 1 for X t as well, since the expected exit time 



cannot decrease, and the Green functions for X t can only grow by a 2d factor since 
the escape probabilites for X t and Y t are at most a 2c? factor apart. 
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Fix x G T, set Tq — 0, and define for each i£N the r.v.'s 

T i =m{{t>T i _ 1 :Y t = x}. 

Let i* = inf {i : T, — oo}. For 1 < i < i* let p t = Ti — Tj_i. We show there are 
positive constants k±, &2 dependent only on d such that 

(2.2) > fcir 2 /logr] > ^. 

For some r > 1, let dB^(x) = {v E V : dr(v,x) = r}. By electrical network inter- 
pretation (see e.g. [DS84J), the probability for a walk beginning at x to hit dB^(x) 
before returning to x is (2d)~ 1 C e f /(r), where C e //(r) is the effective conductance 
from a; to dB^(x). Since T is infinite and connected, for any r there is a connected 
path of r edges from x to dB^(x). By the monotonicity principle, G e //(r) is at 
least the conductance on this path, which is r _1 . Thus the probability to hit some 
y G dB^(x) before returning to x is at least (^dr) 1 . 

Next, let y G dB^(x). By the Carne-Varopoulos upper bound (see |Var85| L 

P [Y t = x\Y = y}< M 1 / 2 exp 

and thus, for some k\(d),k2{d) > and all r > 1, the probability that a walk 
starting at y G dB^(x) does not hit x in the next [kir 2 / logr\ steps, by union 
bound, is greater than &2. Together with our lower bound on the probability that 



we arrive at such a y G dB l r (x), we get (2.2 1 



Next, let g = inf {i : Ti > tz} be the number of visits of Y t to x before hitting 
Z, including t = . g is a geometric random variable with mean G = Gz- Let 



a = h In (4/3) and note that since there is a constant in (2.1l and tz > 1, we can 
assume G > 2. Thus 

P[.9>«G] > (l-G- 1 )" 6 

> (1-G- 1 ) 2 "^ 

> e- 2Q =3/4 

We further assume G > | V || so that aG - 1 > aG/2 and Gf§ > 1. Let A be the 
event that there is an 1 < i < (aG — 1) A i* such that p % > fci (j§G) 2 /log (j§G) 



Note that i* < aG— 1 implies A. Thus by (2.2 1 and the independence of consecutive 
excursions from x, 

( 16\ qG/2 
P [A c ] < (l - < e- 8Q 

which is smaller than 1/4. 
Thus 

P[{g>aG},A]> \. 

This implies the lemma since tz > X)f=i Pii anc ^ f° r & = 53i=i Pi > 

fcG 2 /logG. ' 1 " □ 

Next, we state the fact that a Harnack inequality implies an oscillation inequality. 
For a set of vertices U and a function u, let 



oscu(u) = max;/(tt) — min[/(w). 
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Proposition 2.2. Let x G T and assume that for some r > we have that any 
function h that is non-negative and harmonic on B2 r {x) satisfies (1.5! on B r {x). 
Then for any h that is harmonic on B2 r (x), we have 

1 



(2.3) 



OSC Bt ( x) 



h < 



ch 



-OSC B2r (x)l 



CH + l 

Proof. We quote a proof from chapter 9 of |Tel06j . 

Set v = h — mia. B 2r .( x ) h. Since v is non-negative and harmonic on B2 r {x), the 
Harnack inequality (1.5) is satisfied here, so we have 



whence 



and 



max h 

B r {x) 



max v < ch min v, 

B T (x) B T (x) 



min h < ch [ min h 



B 2r (x) 



B r (x) 



min h 

B 2r (x) 



oscr ( x )h < (ch — 1) min h — min h 

\B r (x) B 2r (x) 



Similarly, we have 



osc Br(x) h < (c H - 1) 



max h 

B 2r (x) 



max h 

B r (x) 



Summing up these two inequalities, we obtain 

1 



osc Br(x) h < - (c H - 1) (osc B2r{x) h ~ osc Br{x) h) , 



□ 



whence (2.3) follows. 

Iterating this on the Harnack assumption [4] we get 

Corollary 2.3. For any a > there is an M(a) such that for any r\ > and for 

all R > R v , if x £ Bfj, r > nR and h is a harmonic function on B M ^ r (x) then 



osc Br{x) h < aosc Bjl 



,h. 



We use this to show regularity of the green function and for expected exit time. 
Lemma 2.4. There is a cc(e) such that for all large enough R, and any x £ 

B(l- e )R \ B eR 

(2.4) G TR (0,x)<c G R 2 - d . 

Proof. Any two vertices in u, v £ -B(i_ c )t? \ B e R can be joined by a path of n < 

C(d)/e overlapping balls B e / 4 (xi), . . . , B e / i (x n ) that are all subsets of Bn_ e /2)R \ 

B(e/2)R such that x\ = u and x n = v. Since G TR (0,x) is positive and harmonic in 

B(i—e/2)R \ B{e/2)Ri Harnack assumption |4] tells us that c H C ^ < u/v < c H ^. Next, 

note that J2xgb r G r R (0,x) = E [tr]. Thus, letting M(R) = max {G TR (0, x) : x € B (1 _ e ) R \ B eR }, 

and using the exit time bound in assumption [3] we get 

\B R \ c H c/€ M < G TR (0,x) < c E R 2 . 

xEBr 



a 
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Lemma 2.5. For any f3 > 0, there is a S > such that for all large enough R, any 
x, y G -B(i_ £ )# \ B e R satisfying \x — y\ < SR also satisfy 



(2.5) 



G TR (0,x)-G TR (0,y) 



< PR 



1-d 



Proof. G TR (0,x) is positive and harmonic in -B(i_ e /2)i?, \ -B( e /2).R and bounded by 
ccR 2 ~ d - We use corollary 2.3 with a — fic^ 1 and r\ = e/ (2M(a)). This gives us 
the lemma with S = q. □ 

Lemma 2.6. For any ft > 0, there is a S > such that for all large enough R, any 
x,y G -B(i_ c )a satisfying \x — y\ < SR also satisfy: 



(2.6) E x [T R ]-E y [T R ] <(3R 2 . 

Proof. Fix x G -B(i_c)i? and for <5 > determined below, fix some y G B$r(x). 
Define the r.v. n to be the random time it takes a walk starting somewhere inside 
B € r(x) to exit B t n(x), and let T2 = tr — ti, i.e. the additional time it takes the 
walk to exit Br. Note that -E^ [t2\ is harmonic in B tR (x), and that from exit time 
assumption [3] for all large R it is bounded by ceR 2 - We use corollary 2^3 with 
a = (3c E 1 /2 and rj = e/M(a), to get that for S < r\ 



E x [r 2 ] - E y [t 2 



< ac E R 2 = PR 2 /2. 



Take 5 < ( J?— ) An so that, again by exit time bound, for all large enough R, 



1/2 



E x [ti] , E y [ti] < j3R /2. Applying the triangle inequality finishes the proof. □ 

3. Domination of Green function 

Let = r NLJ {°> and let P B {-) denote a probability measure on paths in f] 
star ting at that make X t a "blind" simple random walk as defined in subsection 
P B is pushed forward to a measure on by w R {t) as defined in subsection 
To contrast, we call P(-) the Wiener measure on curves corresponding to the 



1.2 



1.4 



Brownian motion B(t) which is the weak limit of wn(t). Thus, for fixed T, is 
the probability space on which P B [w R {t)) converges to P in distribution. Write 
E B [■] and E [■) for the corresponding expectations. 

3.1. Integral Convergence of expected exit time. Since we assume control of 
convergence to Brownian Motion only from 0, we must describe the expected exit 
time from an arbitrary point in the unit ball as a function of Brownian motion that 
starts at 0. We do this by conditioning the Brownian motion to hit a small box 
containing that point and measuring the additional time needed to exit the unit 
ball. 

We denote the first hitting time of a set Z by tz ■ This hitting times may refer to 
the Brownian motion B(t), scaled linearly interpolated walks wn(t), or the discrete 
random walk X t . Another implicit part of this notation is the starting point of 
the walk or curve. The correct interpretation should be evident from context, and 
will be stated otherwise. Some notation used in this section was introduced in 
subsections 11.21 and 11.41 
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Fix T > 0, < 6 < e, u € bi e and let A = A(T) = {w(t) E C$ T : r de(u) < T}. 
A C C|t is the event that the curve hits a small box around u before time T. 

Henceforth, to avoid the complication where a vertex after scaling is in the 
boundary of a box dg(u), we always take u to have rational coordinates, and the 
side length 9 to be irrational. This will suffice as our scale parameter R is a natural 
number. Secondly we take T to be an integer so that a curve w R (t) hits dg(u) until 
T if and only if it hits a vertex in D R g(Ru) until TR 2 . 

Since we are interested in estimating the behavior of X t and not just its interpo- 
lation, we define A* = A*(R,T) = {u e fl : T Dm{Ru) (X t (u)) < TR 2 }. A* C fl is 
the event a vertex in D R g(Ru) is visited until time TR 2 . Thus for all R e N, and 
for all u€.Q, X t (u) e A* w R {u,t) G A. 

Let t + be the first exit time of the unit ball bi after hitting dg(u). Let 
k : C| T -> K be defined as 

fc(w) = U((r + -T d9(u) )AT). 

Analogously, let be the first exit time of after hitting Dng(Ru), and let 

fc» - i?- 2 • 1 A , ((r+ - r DR9(iJu) ) A TR 2 ). 

k(u) is bounded by T and is discontinuous on dA, a set of Wiener measure zero. 
Therefore, by the Portmanteau theorem and our assumption of weak convergence, 

E B [k (w R (t))} ->E[k (B(t))} as R -> oo. 

By the strong Markov property for Brownian motion, we may average over the 
hitting point. 



P(A)- 1 ■ E [k (B(t))} =E[k (B(t)) | A] = E 



E 



B(T dfl(u) ) 



AT] 



where r = Tbj, the first exit time from the unit ball (we start measuring time at 
^( T d e (u)))- Note that R 2 k(wn(u,t)) measures the time that the unsealed interpo- 
lated walk w\(u,t) takes to get from the boundary of dRe(Ru) to b^,, but what 
interests us is the span between the first time that X t (u) takes a value in Dng(Ru) 
to the first time it takes a value in the complement of B R . R 2 k* (X t (u)) measures 
this time. By the strong Markov property for random walks 



P B (A*)- 1 -E B [k* (X t )] = E B [k* (X t ) | A*] 



R~ 2 E B 



E 



y ( T D S( (Bu)) 



[tr A TR 2 } 



If the unsealed interpolated curve w\(u, t) crosses the boundary of djig(Ru), it will 
hit a vertex in D R g(Ru) in less than one time unit. The same is true for exiting 
b R . Thus for all u in our probability space \R 2 k (wr(u, t)) — R 2 k* (X t (u))\ < 2 
and \E B [k(w B {u,t))]-E B [k*(u)]\ < R \. 

By weak convergence, for any fixed T, P B (w R (u,t) e A(T)) -> P(A(T)), and 
since A* <=> A, 



(3.1) R 



-2 rpB 
^0 



E B 



[tr A TR 2 



En 



E, 



w(T d9 (u)) 



[t A T] 



as R 



oo. 



In summary, we have some average on the boundary vertices of D R g(Ru) of the 
function R~ 2 E b (t r A TR 2 ), that is close as we like to an average on dg(u) of a 
Brownian motion's expected time to exit the unit ball. 
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3.2. Integral convergence of Green function. For a fixed T > 0, > and 
u e b 1 let h : Cf, — > K be defined for to(t) € as follows: 

r 

h(w(t))= J l{ w (t)ed B (u),t<T}dt. 
o 

/i(w) measures a curve's occupation time of dg(u) before leaving bi and until time T. 
Since h(B(u>,t)) is bounded by T and is discontinuous on curves whose occupation 
time of <9dg(u) before exit of bi is positive - a set of Wiener measure zero. Thus 
the Portmanteau theorem gives: 

E B [h (w R (t))] -+E[h [B(t))\ as R oo. 

Note that E [h (B(t))] = J g TA T(0, x)dx where g T /\T is the Green function of B(t) 

d„(u) 

killed on leaving the unit ball or when time T is reached. 

Again, we are measuring the time that a linearly interpolated curve spends in 
a set, while we would like to have control over the time the random walk itself 
spends in the set. However, E [h (B(t))} is a continuous function of 6, and for any 
5 > the time the discrete walk spends in D R e(Ru) is eventually sandwiched 
between the time the unsealed interpolated walk wi(u),t) spends in d^g+a) (i?u) 
and dj{(0_a)(i2u). Thus: 

(3.2) R- 2 Yj G ? r at&(°> x )^ J 9rAT(0,x)dxasR^oo, 
xeD m (B.u) de(u) 

where G^ rATR2 is the Green function of the walk from 0, killed on leaving Br or 
when time TR 2 is reached. 

3.3. Pointwise domination of Green function. 

Lemma 3.1. Ve > 3R and r\ > s.t. Vi? > R the following holds: 
xeB (1 _ e)R \B eR =^ \B R \G TR {0,x)> (l + n)E x {T R ). 

This is the main result needed for IDLA lower bound, and will proved be in this 
section. 

First, it is known (see, e.g., [LBG92J p. 2121) for Brownian motion starting at 
zero, and killed on exiting the unit ball, that the Green function g T (0,x) and ex- 
pected hitting time E x (t) are continuous functions with the property that |bi|<jv(0, x) — 
E x (t) descends strictly monotonically to zero as x goes from to 1. This is true 
for any Brownian motion, in particular, B(t), the weak limit of X t on the graph 
starting at 0. Thus for any e > the minimum of the difference between the two 
when \\x\\ € [e, 1 — e] is bounded away from zero. The Lebesgue monotone con- 
vergence theorem implies that g TA% r(0,x) /• g T (0,x) and E x (t A T) /• E x (t) as 
T — > oo. Since all functions involved are continuous and converge monotonically on 
a compact set, by Dini's theorem, the convergence is uniform. Thus we have the 
following uniform bounding of the difference away from zero: 

37(c) >0, Ts.t. Vu,||u|| G [e/2,l-e/2], E u (t A T) + 7 < |b x | 5TAT (0, u). 

Since E u (t A T), <7 t at(0, u) converge uniformly with T, they are uniformly 
equicontinuous in the variable u in the closed interval [e/2, 1 — e/2]. We may 
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then choose a 8 > such that for all large enough T, any average of <?tat(0, •) in 
a box of side 8 with center u. is close to <? t at(0, u) for any u € [e, 1 — e]. We have 
the analogous claim for E u (t A T). Thus: 

Lemma 3.2. For any positive e, there is 7(e) > such that for all large enough T 
and all small enough 8 

Vu,||u|| e [e,l-e], J E v (TAT)d l i(w) + 3 1 <\b 1 \8- d J 5rAT (0,v)dv. 



9d e (u) 



d e (u) 



In the above, d[i is an arbitrary probability measure (total mass one) on the 
boundary of dg(u), while the integral on the right is by rf-dimensional Lebesgue 
measure. 



We apply lemmas 2.5 and 2.6 to get a S > for which (2.6 1 and (2.51 hold with 



/3 = 7 for all large enough R. 7 > is some multiple of 7 from lemma |3.2| that we 
determine later. Set 8 to be small enough so that dg is covered by a ball of radius 



S. Increase T further if necessary so that (1.4) holds with (3 = ?. Fix T and 9 for 



the remainder of the proof. 

We cover b(i_ e ) \ b e by a finite number of open 0-boxes, and so to prove lemma 
3.1 it suffices to prove the implication in the restricted setting of x € Dng(Ru) 
where u is the center of an arbitrary box in our #-net. 

Now, we show the Green function at every point in this box is close to the 
continuous one. 

3.3.1. Pointwise Green function estimate. Let = G Tr ^tr 2 (0, x) and 

xeD R e(Ru) 

let g$ — J g TA T(0,x)dx. By (3.2 1, we have for large enough R: 

d„(u) 



R- 2 Gl - gS 



<10 d \bi\ 



Let a^i — prcatiL)! • ^ density assumption ([2]), for all R large enough 



OtR. 



e-*\\>i\ 



we get 



< 7. Using the triangle inequality on 



R- 2 Gl{a R + 8-^1 ~ a a ) - Mbi| ff 



3-d 



-d\ 



< 7, 



R- 2 Gla R -e-^lgf 



<7 + i?^ 2 G| 



an 



- d |bi| <7(1 + C5^), 

where for the right inequality, we used (2.4 1 to bound G Tr/ \tr 2 (0, %) and (R8) d as 



a bound on the number of vertices. 

We set 8 such that ( |2.5[ ) bounds the difference between the maximum and min- 
imum of G TR (0,x) in D R g(Ru) by jR 2 ~ d for all large enough R. Thus, since 
G TR (0,x) > G TrA tr 2 (0,x) we have for any x 6 D R g(Ru) 



(r- 2 G Tr {Q,x)\B r \ - <r d |bi| 9 /) > - (7(1 + c 5 d ) + R- 2 \B R m 



2-d\ 
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Recall 7 from lemma [372] We now determine 7 so that for all large enough R, for 
any x £ D R g(Ru) 



( 



(3.3) 



R- 2 G TR {0,x)\B R \-6- d \h 1 \ J g TAT (0,v)<h 



\ 



> -7- 



\ d e (u) / 

Next we show that the expected hitting time from any point is close to the contin- 
uous one. 

3.3.2. Pointwise expected hitting time estimate. Recall we chose T such that for 
any x € B R 

\E x [t r ]-E x [t r ATR 2 ] \ < ^R 2 . 
Combining this with ( |3.1[ ) we have for large enough R 

R~ 2 E ^E x (r R )- J E v (t A T)rf^i(v) 



x£D Re (Ru) 



ade(u) 



< 



where for each R, X x are non-negative and sum to one and dfj, is some probability 
measure on ddg(u). 



We set 6 such that (2.6 1 holds with j3 — 7 < 7/2, so for large enough R we have 
for any x £ D R g(Ru) 



R- 2 E x (t r )~ J £ v (TAT)d/x(v) 



dde(u) 



< 7- 



Putting the above together with (3.3 1 and 3.2 for all large enough R, we get that 
for any x £ D R e(Ru) 

G TR (0,x)\B R \ -E x (t r ) > 1 R 2 . 

The above along with our upper bound on E x (t r ) from assumption [3j implies 
lemma 13.11 

4. Lower Bound 
We begin by formally defining the IDLA process. 

Let be a a sequence of independent random walks starting at 0. The 

aggregate begins empty i.e. 1(0) = 0, and X(n) — X{n — 1) U X™, where t' — 
min f {X[ l ^ X(n — 1)}. Thus we have an aggregate growing by one vertex in each 
time step. 

As in |LBG92| . we fix z £ B( 1 _ e ) R \ B eR , and look at the first \B R \ walks. Let 
A = A(z,R) be the event z £ I(\B R \). We show the probability this does not 
happen decreases exponentially with R. 

Let M = M(z,R) be the number of walks out of the first \B R \ that hit z before 
exiting B R . Let L = L(z,R) be the number of walks out of the first \B R \ that hit 
z before exiting B R , but after leaving the aggregate. Then for any a, 

P(A C ) < P(M = L)< P(M < a) + P(L > a). 

We choose a later to minimize the terms. In order to bound the above expression, 
we calculate the average of M and L, 
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E[M} = \B R \P { Tz <T R ). 

E [L] is hard to determine, but each walk that contributes to L, can be tied to the 
unique point at which it exits the aggregate. Thus, by the Markov property, if we 
start a random walk from each vertex in Br and let L be those walks that hit z 
before exiting Br, P(L > n) < P(L > n), and so it suffices to bound P(L > a). 
E[L] is a sum of independent indicators: 

E[L] = p *(tz < t r ). 
xeB R 

Since both M and L are sums of independent variables, we expect them to be 
close to their mean. Our aim now becomes showing that for some 5 > and all 
large enough R. 

(4.1) E[M}> (1 + 26)E[L}. 

By standard Markov chain theory, 

P x (t z < tr) - 

^tr\Z, Z) 

Using this and symmetry of the Green function, it is enough to show for all large 
enough R 

\B R \ G tr (0, z)>(1 + 28) G TR (x, z) = (l + 2S)E Z [r R ] , 

x£B R 

which is lemma I3TT1 

Choosing a = (1 + S)E[L] we write 



(4.2) P(L > (1 + 5)E[L]) < P{L~E[L] > 8E[L] 1/2 a L ) 

(4.3) < 2cxp (-E[L]5 2 /4j . 

In the first line we use that L is a sum of independent indicators, and the variance of 
such a sum is smaller than the mean. The second line is an application of Chernoff 's 



inequality. Similarly, using (4.1 1 



(4.4) 
(4.5) 



P(M < (1 + 6)E[L]) < P(M < 



1 



1 + 



E[M] 



< P(\M — E [M]\ > -E [M]) 



(4.6) < 2cxp(-£;[M](57l6) . 

To lower bound E[L] — E z [tr] (G Tr (z, z))^ 1 we use ( |2.1[ ) to write that for some 
M, 

G TR (z,z) < (ME z [tr\ log E z [t r ]) 1/2 
and since E z [tr] > R we have 

cR 1 / 2 



E 



L(z,R) 



> 



\ogR 
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Together with (4.3 1 and (4.6 1 and summing over all z 6 -B(i- e ),R \ B^r, we get that 
the probability one of the vertices in B( 1 _ e )/j \ B e R is not in Z(|J3r|) is bounded by 
CR d exp(-cR^ 2 /logR). Since the expression is summable by R, by Borel Cantelli. 
this happens only a finite number of times with probability one. So if R' is the 
largest radius for which some vertex in Bfl'(i_ e ) \ -Br/ c is not covered after \Br>\ 
steps, then we possibly have a finite sized hole in the aggregate which will almost 
surely fill up after another finite number of steps. 
Thus we have proved the main theorem |1.1| 

Acknowledgement. Thanks to Itai Benjamini for suggesting the problem, to Gady 
Kozma for many helpful discussions, and to Greg Lawler for suggesting an elegant 
approach. 
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